Collectivity-induced quenching of signatures for shell closures 



00 

o 
o 

(N 

< 
(N 



> 

00 

o 

cn 
od 

o 

oo 
O 



M. Bender, G. F. Bertsch,^ and P.-H. Heenen'''^ 

^ Universite Bordeaux, Centre d'Etudes Nucleaires de Bordeaux Gradignan, UMR5797, F-33175 Gradignan, France 
^CNRS/IN2P3, Centre d'Etudes Nucleaires de Bordeaux Gradignan, UMR5797, F-33175 Gradignan, France 
''Department of Physics and Institute for Nuclear Theory, 
Box 351560, University of Washington, Seattle, WA 98195 
'^Service de Physique Nucleaire Theorique, Universite Libre de Bruxelles, CP. 229, B-1050 Bruxelles, 

''Physics Division, Argonne National Laboratory, 
9700 S. Cass Avenue, Argonne, IL 60439, U.S.A. 
(Dated: 27 August 2008) 

Mass differences are an often used as signature and measure for sliell closure. Using the angular- 
momentum projected Generator Coordinate Method and the Skyrme interaction SLy4, we analyze 
the modification of mass differences due to static deformation and dynamic fluctuations around the 
mean-field ground state. 

PACS numbers: 21.10.Dr, 21.60.Jz, 



I. INTRODUCTION 

How the magic numbers that are well established along 
the line of stability evolve when moving away towards nu- 
clei with large neutron or proton excess is still an unset- 
tled question. Major recent experimental progress now 
allows one to follow magic numbers over long isotopic or 
isotonic chains far from stability. The data for many ob- 
servables point to the disappearance of well established 
magic numbers and the appearance of new ones in nuclei 
which sometimes are still close to stability, and far from 
the drip lines. 

From a theoretical point of view, the situation is still 
unclear. The interactions used in shell model calculations 
have had to be revised in order to be able to reproduce 
the new experimental data and to describe the evolution 
of shells, especially in nuclei in the pf shell P, [2, S 0i S 
i. 

An analysis in terms of the self-consistent mean-field 
method performed well before most experimental data 
became available 0] has shown that, when approach- 
ing the drip lines, the increasing diffuseness of the nu- 
clear surface for neutrons and the interaction between 
bound orbitals and the continuum affect the shape of 
the mean- field potentials of protons and neutrons. As 
a consequence, shell structure is modified, known gaps 
are quenched, and new ones might open. However, the 
coupling between bound orbitals and the continuum can 
be expected to be an important mechanism for nuclei in 
the immediate vicinity of the drip line only, when the 
Fermi energy is approaching zero. With the exception 
of the case of the lightest nuclei with mass below 30, 
there are no experimental data for such nuclei, and it is 
not obvious that the coupling with the continuum is an 
effect which has an impact on the known anomalies in 
shell structure. By contrast, the diffuseness of the neu- 
tron density might play a role in unstable nuclei, mainly 
through a reduction of the spin-orbit potentials for both 
protons and neutrons [1,111 . However, this effect can be 
expected to appear gradually, and cannot not be respon- 



sible for those anomalies in shell structure where nuclear 
structure is suddenly changing when adding or removing 
a few nucleons only. 

There are several observables that are used as possible 
signatures to put the evolution of shell closures in evi- 
dence when following isotopic or isotonic chains. From 
a theoretical point of view, the simplest one is the two- 
nucleon separation energy in even-even nuclei, which can 
be easily calculated from the total binding energy. A 
jump in the two-nucleon separation energy is a direct in- 
dication of a sudden increase in the ground-state binding 
energy of a given nucleus, and often has a shell closure 
as its origin. Other quantities that are frequently ana- 
lyzed and which arc not directly related to the ground- 
state properties, but to the gaps that separate the differ- 
ent shells are provided by excitation spectra in even and 
odd-mass nuclei. At least in the context of mean-field- 
based models these are more difficult to calculate 
and have larger associated theoretical uncertainties. Still, 
any conclusion about the evolution of shell structure re- 
quires clear clear trends from a variety of observables that 
probe complementary aspects. In particular, it should be 
supported also by data on electromagnetic moments and 
transition probabilities which provide information on the 
structure of the ground state and the excited-state wave 
functions. 

Although there is no formal justification for this prac- 
tice, the discussion of shell evolution is often conducted 
by theorists and experimentalists alike as if the following 
three quantities were equivalent: 

1 . single-particle energies in a spherical self-consistent 
mean field. They are then the eigenvalues of a 
single-particle Hamiltonian and the corresponding 
wave function is a single Slater determinant, i.e. a 
pure state; 

2. effective single-particle energies as defined within 
the interacting shell model. They should in princi- 
ple result from an average over a huge number of 
configurations 0, [H, [II Q E [H, [l3 ; 
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3. two- nucleoli separation energies. 

The last quantity is the only one which can be directly 
related to experimental data. Let us recall its definition: 

S2p{Z,N) = E{Z ~2,N)^E{Z,N) (1) 
S2n{Z,N) = E(Z,N -2)^E{Z,N), (2) 

where E{Z^ N) is the total binding energy of the nucleus 
Z,N. 

In the highly idealized pure HF case where the Koop- 
man theorem is valid [H, [l^ , the two-nucleon separation 
energy is equal to —2 times the energy of the doubly- 
degenerate orbital occupied by the two additional nucle- 
ons. This is never exactly the case in nuclear physics. 
First, pairing effects require a change in the energy of 
the Fermi level, resulting through self consistency to a 
change of all orbits. Even without pairing correlations, 
the two extra particles have an effect on the core and may 
induce at least a small rearrangement of all the orbits, or 
in the worst case, the appearance of deformations. 

It is the purpose of this paper to analyze the role of 
static and dynamical quadrupole correlations on the sys- 
tematics of S2q{Z, N) in the context of mean-field based 
methods. Our analysis relies on the results of Ref. [20l.[2]| 
where systematic calculations of the masses of even-even 
nuclei were presented. The method allows one to intro- 
duce the correlations due to symmetry restorations (par- 
ticle number and angular momentum) performed on axial 
mean- field configurations and to add a mixing of the pro- 
jected configurations with respect to the axial quadrupole 
moment in the framework of the discretized Generator 
Coordinate Method (GCM). Since the starting point of 
the method is a mean-field calculation, the results give 
access at the same time to the spherical single-particle 
spectra of all even-even nuclei and to the two-nucleon 
separation energy of states including beyond mean-field 
correlations. 

In the next section, we recall the main features of our 
method. We then discuss in details the evolution of the 
two-nucleon separation energy along Z or = 50. Fi- 
nally, we compare the results obtained for all isotopic 
and isotonic chains with the three different steps of our 
method to the experimental data. 

II. THE MODEL 

The method used to calculate binding energies for the 
ground states of even nuclei is described in detail in 
Rcfs. [2^ [21} (referred to as paper I in the following) . 
In this analysis, we use the energies as tabulated in [22l |. 

To summarize the key features of our method, its start- 
ing point is a set of mean-field calculations including a 
constraint on the axial quadrupole moment. As effec- 
tive interaction we employ a Skyrme energy density func- 
tional, the SLy4 parametrization [23| . for the mean-field 
channel, and a density-dependent, zero-range interaction 
in the pairing channel. Two sets of correlations beyond 



the mean-field are introduced. First, the deformed wave 
functions are projected on both fixed particle numbers 
and on angular momentum. This projection corresponds 
to mixing degenerate mean-field wave functions which 
differ by a spatial rotation and to introduce fluctuations 
around the orientation of the mean-field. In a collective 
model terminology, these correlations would be called ro- 
tational correlations. A second step of our method con- 
sists in the mixing of projected mean-field wave functions 
as a function of their axial quadrupole moment. This 
corresponds to a vibrational correction in the collective 
model language. Our final wave function has the form: 

\JM,.)^Y.fj^'^PMoPNPz\q). (3) 

The ket \q) is a (paired) mean-field state of axial 
quadrupole deformation q. The particle numbers and an- 
gular momentum quantum number are restored thanks to 
the operators Pjv and Pz which project on good neutron 
and proton numbers, and P^jq which projects on angular- 
momentum J with the z component M in the laboratory 
frame The projected wave functions are labeled by 
the index q which indicates from which mean-field state 
they are obtained, although an intrinsic quadrupole mo- 
ment has no direct physical meaning for a wave function 
with good angular momentum and can only be defined 
with the help of the collective model, in the same way as 
it is from experimental data. The weights fj^^{q) defining 
the mixing of the projected wave functions with respect 
to q are obtained by variation of the total energy. 

The mean-field wave functions are generated by treat- 
ing the pairing correlations using the Lipkin-Nogami 
(LN) prescription. This enforces the presence of pair- 
ing in all nuclei and for all quadrupole moments, even at 
magic numbers, which is necessary to ensure the continu- 
ity of the mean-field wave functions required to perform 
the configuration mixing. Binding energies, however, are 
always recalculated after a projection on particle num- 
ber 

We stress that there are no assumptions made in the 
model about the amplitude of the quadrupole fluctua- 
tions introduced into the calculations. Depending on the 
structure of a nucleus, this amplitude either corresponds 
to a small vibration around a pronounced minimum, to a 
large-amplitude motion in a soft and wide potential well, 
or to the mixing of several states around coexisting min- 
ima in the deformation energy surface. Let us recall also 
that our method is fully variational and that the energy 
of the ground state is lowered at each successive stage of 
the calculation. 

Finally, one must note that the Skyrme interaction is 
not a force, but an energy density functional. It requires 
some care to be used beyond a mean-field approach, and 
there are still some open questions concerning this long- 
standing practice. The density-dependent terms are gen- 
eralized using the standard prescription that the den- 
sity is replaced by the transition density in the density- 
dependent terms. This is the only prescription that 
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guarantees various consistency requirements of the en- 
ergy functional pBI . [26| . However, this procedure may 
lead to problems that have been put into evidence re- 
cently [23,[H,[1^. The numerical procedure that we have 
used (see Paper I for details) appears to be safe in this 
respect and no dramatic signs of problems have shown 
up when varying the number of discretization points in 
deformation. In any case, since we arc interested here in 
trends as a function of N and Z , errors of the order of 
200 to 300 keV would not affect our conclusions. 

We will now compare results obtained from three wave 
functions that successively add quadrupolc correlations: 

1. spherical mean field states \q = 0); 

2. the mean- field minimum in the space of axial 



reflection-symmetric deformations 
might be spherical; 



which 



3. the ground state obtained after configuration mix- 
ing by the generator coordinate method of J = 
projected axial quadrupole. We refer to these wave 
functions in the following as projected GCM. 



> 

CO 





-2 
-4 
-6 
-8 
-10 
-12 
-14 
-16 



-10 



y -11 



o 



a. 



-12 
-13 



-14 



-15 
-16 




O spherical mean-field 
□ deformed mean-field 
• J=0 GCM 



experiment 



Z=52 




oooQooOaDD 

Z=50 



50 52 54 56 58 60 62 64 66 68 70 72 74 76 78 80 82 84 

Neutron Number N 



III. RESULTS 

We will first illustrate the key points of the analysis 
of two-nucleon separation energies as measures of shell 
closures for the proton Z ~ 50 and neutron A'^ = 50 
shells. In a second step, we will then analyze the global 
systematics of two-nucleon separation energies for even- 
even nuclei across the nuclear chart. 

Recent experimental progress has provided many 
new direct, high-precision measurements of nuclear 
masses [13, HH, which allow comparison of calculations 
with data along chains that extend to exotic nuclei. 
Experimental masses used here are either taken from 
AME2003 [33], where we exclude extrapolated values, or 
from Ref. j30l . [Slj when error bars are smaller or when 
the AME2003 values from indirect measurements do not 
agree with new directly measured ones. 

A. The Z = 50 isotopic chain 

Figure [T] illustrates the difficulties encountered when 
relating two-nucleon separation energies to the gaps in 
the single-particle spectrum with the example of the 
chain of even-even Sn, Z = 50, and Te, Z = 52, isotopes. 
The lower panel shows the evolution of the eigenvalues 
of the spherical mean- field single-particle Hamiltonian e^. 
Note that all Sn isotopes have a spherical ground state 
at the mean-field level. 

The open blue circles, plotted in both panels, represent 
the two-proton separation energies S2p for iV = 50 and 

= 52 from spherical mean-field calculations. 

The magnitude of the Z = 50 gap in the single-particle 
spectrum is fairly independent of N in this case. The 



FIG. 1: (Color online) Top: Single-particle spectra of the pro- 
tons for the chain of Sn {Z = 50) isotopes and —1/2 times 
the two-proton separation energy for the Z = 50 and Z = 52 
isotopic chains (see text). Bottom: Comparison of the gaps 
between -S2p(Z = 50,iV)/2 and -S2p{Z = 52,N)/2 ob- 
tained at three levels of approximation with the experimental 
data. A linear trend equal to {N - 82) [S2piZ = 50, A = 
50) - S2p{Z = 50, A = 82)]/2 using the spherical mean-field 
values has been subtracted from all —S2p{Z, N)/2, such 
that the distance between the curves for a given is un- 
changed. In particular, the plotted values at A = 82 are 
not afi^ected, and the plotted spherical mean-field values for 
{Z = 50, iV = 50) and {Z = 50, iV = 82) are identical. 



spacing between the — S'2p/2 for Z = 50 and 52 is a poor 
measure of the true shell gap, even when — S'2p/2 values 
and the single-particle levels are calculated in the same 
framework and the same shape is enforced for all nuclei. 
Indeed, the — S'2p/2 spacing turns out to be smaller for 
all A^ values. This result can be partly attributed to 
our treatment of pairing by the LN prescription. This 
method has the effect of moving the single-particle states 
away from the Fermi level. However, the amplitude of 
this change is of the order of 200 to 300 keV at most and 
does not account for the difference between the single 
particle levels and — S'2p/2. 

Still, the evolution with A^ of the spherical mean-field 
—Sipjl values follows closely that of the single-particle 
levels: the difference between their values for Z = 50 and 
Z = 52 is nearly constant. 

In the lower panel of this Fig. [2 the — S'2p/2 values for 
both isotopic chains are shown for spherical, deformed, 
and projected GCM calculations together with the exper- 
imental data. A global descending linear trend with A^ 
has been taken out. The —Sipjl values from spherical 
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FIG. 2: (Color online) Same as Fig. [T] but for neutrons. In 
the lower panel, a global descending trend (A'' — 50) [S2n{Z — 
28, N = 50) - S2n{Z = 50, iV = 5Q)]/2 using the spher- 
ical mean-field 5*2^ values has been subtracted from all 
— S2n{Z,N)/2. The overall reduction of the distance be- 
tween S2n{Z, N = 50) and S2n{Z, N = 52) when approaching 
Z = 28 is clearly visible. 



mean-field calculations differ tlie most from the actual 
data. Tfiey are smaller for Z = 50 and larger for Z — 52 
than tlie data. In particular, the distance between the 
experimental iS'2p values at Z = 50 and at Z = 52 is not 
constant; it decreases when going away from N = 82, 
and starts to grow when approaching N = 50. 

A significant improvement is achieved by allowing the 
mean-field to deform. As outlined in paper I, some of 
the Cd (Z = 48) and Te {Z = 52) isotopes gain up 
to 2 MeV in binding energy thanks to deformation, and 
this increases S2p for Z = 50 and decreases it for Z = 52. 
Despite this effect of deformation, none of the nuclei in- 
volved in the calculation of the S2p can be classified as 
well-deformcd: they exhibit a transitional pattern, with 
soft deformation energy surfaces, and present shallow de- 
formed minima in some cases. This soft topography has a 
significant effect when the fluctuations around the mean- 
field minima are taken into account. The Cd and Te 
isotopes are usually softer than the Sn ones for a given 
A^- value, in such a way that Cd and Te gain more dy- 
namical correlation energy. This increases the S2p values 
for Z ~ 50, but reduces them for Z = 52. 

Specifically, for these two isotopic chains, the change 
in the S2p energies brought about by the dynamical 
quadrupole correlation energy is often large when the 



change due to the static deformation energy is small 
and vice versa. The resulting effect is that the pro- 
jected GCM curves are much smoother than the de- 
formed mean-field ones, and closely follow the experi- 
mental data. The only discontinuity left in Fig. [T] is at 
the neutron shell closure A^ = 82, where the distance 
between S2p{N, Z = 50) and S2p{N, Z = 52) is largest, a 
phenomenon coined as "mutually enhanced magicity" in 
the literature [H, [H, [sl] . From this result one can con- 
clude that, within the present theoretical framework, the 
variation of the gap in the S2p values as a function of A^ is 
not related to a change in the underlying spherical shell 
structure, but to the variation of the energy contribution 
brought by static and dynamic quadrupole correlations 
to the masses of the nuclei that enter the calculation of 
S2p- The same qualitative behavior has been found for 
proton shells at Z = 82 in Refs. 0, [SI, [Hi. 



B. The = 50 isotopic chain 

The same analysis for the chains of A" = 50 and 52 
isotones is performed in Figure [2l Recent mass measure- 
ments [31[ allow one to follow S2n down to Z = 30. Qual- 
itatively, the results arc quite similar to those discussed 
for the Z — 50 and 52 chains above, with two notable ex- 
ceptions. First, the size of the A^ = 50 gap in the single- 
particle spectrum is slowly and continuously decreasing 
from 6.1 MeV in i°"Sn to 4.4 MeV in ^^Ni, which consti- 
tutes "real" shell quenching and is a consequence of the 
increasing effect of the surface on neutron levels that are 
pushed up in the potential well. The spin-orbit splitting 
between the neutron gg /2+ a-nd 177 /2+ levels is slightly in- 
creased in our calculation when going from ^''"Sn down 
to ^''Ge, the lightest isotone where both are bound. The 
striking difference with respect to the protons in the Sn 
chain discussed above can be related to the impact of the 
Coulomb barrier: its presence for protons suppresses the 
diffuseness of the wave functions. Second, several fluctu- 
ations are superimposed on this smooth global trend, one 
of which slightly opens the gap at the proton sub-shell 
closures Z = 38 and Z = 40 and closes it for smaller 
Z values, such that the calculated gap is smallest for 
Z 32, in agreement with experiment [3l] |. when tak- 
ing also the experimental S2n value of the odd nucleus 
(Z = 31, A^ = 52) into account. 

A comparison of predictions of various mean-field and 
other mass models for the A^ = 50 isotonic chain with 
the most recent data can be found in Ref . [311 . 



C. Global systematics 

Figure[3]provides the two-proton separation energy S2p 
for isotonic chains (now without the factor —1/2) in all 
but the lightest even-even nuclei calculated in [23, [2l| . 
Together with the experimental data, the plots give the 
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FIG. 3: (Color online) Two proton-separation energy for 
even-even nuclei. Lines connect nuclei in isotonic chains. 



results for spherical and deformed mean-field calculations 
and for the projected GCM ground states. 

The energies from spherical mean-field calculations 
(top panel) nearly follow straight lines, with a slope that 
slowly decreases with mass. Three gaps, corresponding 
to the magic numbers Z = 28, 50 and 82, are clearly 
visible. Their widths do not vary significantly with neu- 
tron number. Most importantly, there is no sign of shell 
quenching in the S2p trends, consistent with the gaps in 
the single-particle spectra that remain fairly constant as 
well. 

The inclusion of static mean-field deformations has a 
significant effect on the curves. Compared to the spheri- 
cal results, the magic gaps in the S2p values are reduced 
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FIG. 4: (Color online) Two-neutron separation energy S2n for 
even-even nuclei. Lines connect nuclei in isotopic chains 



in many instances. As in the examples of the Z = 50 
and = 50 gaps, the magic nuclei remain spherical, but 
not all of their neighbors with ±2 protons are. The slope 
of the S2p curves is also modified between the spherical 
shells: the energy gain due to deformation increases 
above magic numbers and decreases it below. Moreover, 
several small local gaps open between the spherical shell 
closures, and indicate deformed shell closures. 

Correlations beyond mean-field amplify the effects 
brought by static deformation: the gaps at the magic 
numbers are reduced further. Confirming the discussion 
for the Z = 50 chain, deformation and dynamical cor- 
relations reduce the gap in the S2p energies at all shell 
closures to about half their value in spherical mean-field 
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calculations. Also, the effect on energy of dynamical cor- 
relations is the largest at shell closures, where it rapidly 
changes from nucleus to nucleus, while it is nearly con- 
stant for open shell nuclei [2l|. It is quite gratifying to 
see that adding the correlations brings the calculated val- 
ues close to experiment. Note that a particularly strong 
quenching is visible for Z = 82. The three calculations 
presented here arc based on the same effective interaction 
which predicts a spherical single-particle spectrum with 
gaps much larger than those in the experimental S'2p en- 
ergies that vary only by a few 100 keV with the neutron 
number. Hence, the quenching observed in the experi- 
mental S^p values results from the effects of deformation 
and of beyond mean-field correlations on the total bind- 
ing energy and could be called "collectivity enhanced" 
quenching of the two-nuclcon separation energies at shell 
closures, rather than quenching of the gaps in the spher- 
ical single-particle spectra. 

Figure |4] shows the two-neutron separation energies 
Sin for isotopic chains. Qualitatively, one finds a sim- 
ilar behavior as for the S^p curves, and the quenching of 
the 5*2,1 values from quadrupole correlations is even more 
pronounced. However, the overall agreement with exper- 
iment of the Sip curves from the J = projected GCM 
is better than the one for 5*271, as was already pointed 
out and discussed in detail in [23, [2l| . The Skyrmc in- 
teraction SLy4 used here (as all others we tested) seems 
to systematically overestimate neutron shell gaps, while 
proton shell gaps are described better. 

The inclusion of beyond-mean field correlations has a 
smoothing effect on the variation of the two-particle sep- 
aration energies as a function of or Z. The often un- 
realistic local fluctuations obtained without these corre- 
lations in our calculations or in those of deformed mean- 
field calculations by other groups [37| , or from mean- field- 
based mass models such as the microscopic-macroscopic 
method [11] or the Brussels HFB mass fits [13 are to 
a large extent suppressed. Thus, the correlation energy 
added by fluctuations of the wave functions around the 
mean-field minima compensates the too abrupt shape 
changes (spherical to deformed, or prolate to oblate) ob- 
tained within a pure mean-field approach. 



IV. SUMMARY 

The structural changes that occur in nuclei when go- 
ing along isotopic or isotonic chains have been discussed 
on the basis of two-nuclcon separation energies. We have 
shown that the effect that has been referred to as shell 
quenching and has sometimes been attributed to a re- 
duction of the spherical gaps far from stability is ade- 
quately described by the introduction of dynamic col- 
lective quadrupole correlations. Much smoother trends 
with fewer kinks and discontinuities are obtained than 
in deformed mean- field calculations. In many instances, 
our calculated values are close to the data. Looking 
at the shell closures where the Siq values exhibit dis- 



continuities, static mean-field deformations and dynami- 
cal correlations decrease systematically the amplitude of 
these gaps, and reduce them far from stability. Both 
effects are not related to a reduction of the spherical 
shell structure, rather both underline the importance of 
fluctuations around single mean-field configurations for 
a high-precision description of nuclear masses. Studies 
performed in restricted mass regions JO, 21,, ,36,, 411, 4lj 
have arrived at similar conclusions. 

This study addresses also two present-day questions 
concerning mean-field models and effective interactions 
(or density functionals). For a long time, single particle- 
energies have been viewed as being not directly con- 
nected with experimental data because they are strongly 
renormalized by the coupling to different kinds of vibra- 
tions [Tsl. [4^ . 0, HI]. However, it has been proposed 
recently [i^ that this statement is too strong and that 
single-particle energies should be used in the fitting pro- 
cedure of energy density functionals. This point has been 
convincingly demonstrated using specific data. Our anal- 
ysis shows that one must carefully select the data relevant 
for adjusting single-particle levels and that the quenching 
of 2-particlc separation energies should not be related to 
single particle properties. Our study also raises the ques- 
tion of whether all correlations should be included in an 
energy-density functional, or if some of them should be 
treated explicitly, beyond a mean-field approach. At the 
minimum, our study shows the interest of treating ex- 
plicitly the correlations associated with fiuctuations in 
collective degrees of freedom. 

Finally, let us note that the relation between single- 
particle energies and the evolution of 2-particle separa- 
tion energies along A" or Z is not a simpler problem in the 
shell model context. The definition which corresponds 
the most closely to a spherical mean-field single particle 
energy is the effective single-particle energy of Caurier et 
al. [ll|, [l^ , which includes the monopole shift of single- 
particle energies, i.e. the terms of the interaction that 
correspond to the spherical mean field. A definition more 
directly related to nucleon transfer reactions is based on 
ccntroids of the strength distribution [3, [H, [IB, Il7| . 
Both are obtained from the averaging over correlated 
states which compromises their interpretation in terms of 
a pure configuration, and arc neither a directly observ- 
able quantity. The main message of this contribution is 
probably that one should not try to carry out too indirect 
comparisons with experimental data and that the notion 
of single-particle states should be limited to the very few 
cases where correlation effects are weak. 
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